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THREE-DIMENSIONAL SIMULATIONS OF INFLOWS IRRADIATED BY A PRECESSING ACCRETION DISK 
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ABSTRACT 

We present three-dimensional (3-D) hydrodynamical simulations of gas flows in the vicinity of an 
active galactic nucleus (AGN) powered by a precessing accretion disk. We consider the effects of the 
radiation force from such a disk on its environment on a relatively large scale (up to ~ 10 pc). We 
implicitly include the precessing disk by forcing the disk radiation field to precess around a symmetry 
axis with a given period (P) and a tilt angle (8). We study time evolution of the flows irradiated by 
the disk, and investigate basic dependencies of the flow morphology, mass flux, angular momentum on 
different combinations of O and P. As this is our first attempt to model such 3-D gas flows, we consider 
a simplest form of radiation force i.e., force due to electron scattering, and neglect the forces due to 
line and dust scattering/absorption. Further, the gas is assumed to be nearly isothermal. We find the 
gas flow settles into a configuration with two components, (1) an equatorial inflow and (2) a bipolar 
inflow/outflow with the outflow leaving the system along the poles (the directions of disk normals). 
However, the flow does not always reach a steady state. We find that the maximum outflow velocity 
and the kinetic outflow power at the outer boundary can be reduced significantly with increasing 0. 
We also find that of the mass inflow rate across the inner boundary does not change significantly with 
increasing 0. The amount of the density- weighted mean specific angular momentum deposited to the 
environment by the precessing disk increases as P approaches to the gas free-fall time (iff), and then 
decreases as P becomes much larger than tg. Generally, the characteristics of the flows are closely 
related to a combination of P and 0, but not to P and individually. Our models exhibit helical 
structures in the weakly coUimated outflows. Although on different scales, the model reproduces the 
Z- or S- shaped density morphology of gas outflows which are often seen in radio observations of 
AGNs. 

Subject headings: accretion, accretion - disks ~ galaxies: jets - galaxies: kinematics and dynamics- 
methods: numerical - hydrodynamics 



1. INTRODUCTION 

Powered by accretion of matter onto a super massive 
(10^-101° Mq) black hole (SMBH), Active Galactic Nu- 
clei (AGN) release large amount of energy (e.g., Lynden- 
Bell 1969) as electromagnetic radiation (10^''-10^'*Lq) 
over a wide range of wavelengths, from the X-ray to the 
radio. The very central location of AGN in their host 
galaxies indicates that the radiation from AGN can play 
an important role in determining the physical character- 
istics (e.g., the ionization structure, the gas dynamics, 
and the density distribution) of their surrounding envi- 
ronment in different scales i.e., from the scale of AGN it- 
self to a lager galactic scale, and even to an inter-galactic 
scale (e.g., Quilis et al. 2001; Dalla Vecchia et al. 2004; 
McNamara et al. 2005; Zanni et al. 2005; Fabian et al. 
2006; Vernaleo & Reynolds 2006). The feedback process 
of AGN in the form of mass or energy outflows, in turn, 
is one of key elements in galaxy formation/evolutionary 
models (e.g., Ciotti & Ostriker 1997, 2001, 2007; Silk & 
Rees 1998; King 2003; Begelman & Nath 2005; Hopkins 
et al. 2005; Murray et al. 2005; Sazonov et al. 2005; Silk 
2005; Springel et al. 2005; Brighenti & Mathews 2006; 
Fabian et al. 2006; Fontanot et al. 2006; Thacker et al. 
2006; Wang et al. 2006, Tremonti et al. 2007). 

Although the AGN outflows can be driven by magne- 
tocentrifugal force (e.g., Blandford & Payne 1982; Em- 
mering et al. 1992; Konigl & Kartje 1994; Bottorff et al. 
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1997) and thermal pressure (e.g., Weymann et al. 1982; 
Begelman et al. 1991; Everett & Murray 2007), it is the 
radiation force from the luminous accretion disk that is 
most likely the dominant force driving winds capable of 
explaining the blucshifted absorption line features often 
seen in the UV and optical spectra of AGN (e.g., Shlos- 
man et al. 1985; Murray et al. 1995; Proga et al. 2000; 
Proga & Kallman 2004). In reality, these three forces 
may interplay and contribute to the dynamics of the out- 
flows in AGN in somewhat different degrees. 

Another complication in the outflow gas dynamics is 
the presence of dust. The radiation pressure on dust can 
drive dust outflows, and their dynamics is likely to be 
coupled with the gas dynamics (e.g., Phinney 1989; Pier 
& Krolik 1992; Emmering et al. 1992; Laor & Draine 
1993; Konigl & Kartje 1994; Murray et al. 2005). The 
AGN environment on relatively large scales (10^ — 10"^ pc) 
is known to be a mixture of gas and dust (e.g. Antonucci 
1984; Miller & Goodrich 1990; Awaki et al. 1991; Blanco 
et al. 1990; Krolik 1999); however, in much smaller scales 
(<^ 10 pc) one does not expect much dust to be present 
because the temperature of the environment is high (> 
10"* K). Concentrating on only the gas component, the 
dynamics of the outflows in smaller scales was studied by 
e.g.. Arav et al. (1994), Proga et al. (2000) in 1-D and 
2-D, respectively. 

Radio observations show that a significant fraction of 
extended extragalactic sources display bending or twist- 
ing jets from their host galaxies. For example, Florido 
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et al. (1990) found that - 11% of their sample (368 ob- 
jects) show anti-symmctrically bending jets (S-shapcd or 
Z-shaped morphology) while ~ 9% show the symmet- 
rically bending jets (U-shaped morphology). Similarly 
Hutchings et al. (1988) studied the morphology of the 
radio lobes from 128 quasars (with z < 1), and found 
that 30% of the sample show a sign of bending jets. The 
bending and misalignment of jets arc also observed in 
parscc scales in compact radio sources (e.g. Linficld 1981; 
Appl et al. 1996; Zensus 1997). Examples of the radio 
maps displaying the S- or Z-shaped morphology of jets 
can be found in e.g. Condon & Mitchell (1984), Hunstead 
et al. (1984), and Trcmblay ct al. (2006). 

Using the data available in literature, Lu & Zhou 
(2005) compiled the list (see their Tab. 1) of 41 known 
extragalactic radio sources which show an evidence of jet 
precession, along with their jet precession periods (P) 
and the half-opening angle {4>) of jet precession cones. 
According to this list, a large fraction (67%) of system 
has rather small half-opening angles, i.e., ip <~ 15°. 
A large scatter in the precession periods are found in 
their sample; however, most of the precession periods are 
found in between lO"' and 10® yr (see also Roos 1988). 
Note that the precession periods are usually too long to 
be determined directly by variability observations. Typi- 
cally the precession periods are found by fitting the radio 
map with a kinematic jet model (e.g., Gower et al. 1982; 
Veilleux et al. 1993). ' Interestingly, Appl et al. (1996) 
showed that a typical precession period of tilted massive 
torus around SMBH is ~ 10*^ yr. 

The S- and Z-shaped morphology seen in the observa- 
tions mentioned above can naturally explained by pre- 
cessing jets. Further, the precessing of jets can occur 
if the underlying accretion disk is tilted (or warped) 
with respect to the symmetry plane. There are at 
least five known mechanisms that can causes warping 
and precession of in accretion disks (1) the Bardeen- 
Petterson effect (Bardeen & Pettcrson 1975; see also 
Schreier et al. 1972; Nelson & Papaloizou 2000; Frag- 
ile & Anninos 2005; King et al. 2005), (2) tidal interac- 
tions in binary BH system (e.g., Roos 1988; Sillanpaa 
et al. 1988; Katz 1997; Romero et al. 2000; Caproni 
et al. 2004), (3) radiation-driven instability (e.g., Pettcr- 
son 1977; Pringle 1996; Maloney et al. 1996; Armitage 
& Pringle 1997), (4) magnetically-driven instability (Aly 
1980; Lai 2003), and (5) Disk-ISM interaction (e.g., 
Quillen & Bower 1999). Using a small sample of AGN, 
Caproni et al. (2006) examined whether mechanisms (1)- 
(4) are capable of explaining the observed precession pe- 
riods. Similarly Tremblay et al. (2006) searched for a 
possible cause of disk precession and warping of the FR I 
radio source 3C 449 using mechanisms (2), (3) and (4) 
above. In general, it is very difficult to determine the ex- 
act cause of disk /jet precession for a given AGN system 
because of large uncertainties in model parameters and 
observed precession periods (which are also often model 
dependent). 

Kochanek & Hawley (1990) presented a hydrodynam- 
ical simulation of jet propagation along the surface of 
an axisymmetric hollow/cone to approximate a jet with 
fast precession; however, intrinsically non-axisymmetric 
nature of the dynamics of jet precession requires the 
problem to be solved/simulated in 3-D. Hydrodynamical 
simulations of extragalactic radio sources with precessing 



jets in full 3-D have been performed by e.g.. Cox et al. 
(1991), Hardee & Clarke (1992), Hardee et al. (1994), 
Typically, in these models, the jets are driven at the ori- 
gin by a small-amplitude precession to break the sym- 
metry and excite helical modes of the Kelvin-Helmholtz 
instability. Careful stability/instability analysis of such 
simulations has been presented by Hardee et al. (1995). 
The effect of magnetic field has been also investigated 
by e.g., Hardee & Clarke (1995) while the effect of opti- 
cally thin radiative cooling on the Kelvin-Helmholtz in- 
stability has been investigated by e.g., Xu et al. (2000). 
Precession of relativistic jets in 3-D with or without mag- 
netic field has been also studied (e.g., Hardee ct al. 2001; 
Hughes et al. 2002; Aloy et al. 2003; Mizuno et al. 2007). 
On much larger scales, Sternberg & Soker (2007) studied 
the effect of precessing massive slow jets onto the inter- 
galactic medium (IGM) in a galaxy cluster, and found 
such jets can inflate a fat bubble in the IGM. In the mod- 
els mentioned above, jets themselves are injected on small 
scales, and the jet propagations are studied. However, it 
is also possible to model a self-consistent production of 
a jet and its subsequent propagation. For example, a 
jet can be produced from an infalling matter by radia- 
tion pressure due to a luminous accretion disk (e.g., see 
Proga 2007; Proga et al. 2007, for axisymmetric cases). 

Regardless of the exact cause of disk/jet precession, 
the observations (e.g. Florido ct al. 1990; Hutchings et al. 
1988) suggest that a significant fraction of AGN contain 
warped or precessing disks. One might expect the details 
of the radiative feedback processes in such systems are 
different from the ones predicted by axi-symmetry mod- 
els (e.g. Proga et al. 2000; Proga 2007; Proga et al. 2007). 
If they differ, then by how much? In this paper, we ex- 
plore the effects of disk precession on the gas dynamics in 
the AGN environment by simulating the outflows driven 
by the radiation force from a luminous precessing accre- 
tion disk around a SMBH. Specifically, we will examine 
how the mass-accretion rate, the outflow powers (kinetic 
and thermal), the morphology of the flows, and the spe- 
cific angular momentum of the gas are affected by the 
presence of a precessing disk and its radiation field. This 
is our first step toward a full extension of the axisym- 
metric radiation-driven wind model of Proga (2007) to a 
full 3-D model. 

In the following section, we describe our method and 
model assumptions, and we give the results of our 3-D 
hydrodynamical simulations in § 3. Our conclusions are 
summarized in § 4. 

2. METHOD 

Our basic model configuration is shown in Figure 1. 

The model geometry and the assumptions of the SMBH 
and the disk are very similar to those in Proga (2007). 
In Figure 1, a SMBH with its mass Mbh and its 
Schwarzschild radius rs = 2GMbh/c^ is placed at the 
center of the cartesian coordinate system (a;, y, z). The 
X-ray emitting corona regions is defined as a sphere with 
its radius shown in the figure. The geometrically- 

thin and optically-thick flat accretion disk (e.g.. Shakura 
& Sunyaev 1973) is placed near the x-y plane. In case 
for an axisymmetric model, the 2;-axis in the figure be- 
comes the symmetry axis, and the accretion disk is on the 
x-y plane. To simulate the disk precession, we assume 
that the angular momentum (Jd) of the accretion disk 
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is tilted from the z-axis by an angle O. In other words, 
the aceretion disk is assumed to be tilted by O from the 
x-y plane. Further, the accretion disk hence its angular 
momentum Jd is assumed to precess around the .z-axis 
with the precession period P. The 3-D hydrodynamic 
simulations will be performed in the spherical coordinate 
system (r, 9, </)), and in between the inner boimdary 
and the outer boundary Tq. The poles of the spherical 
coordinate system coincides with the 2;-axis. The radia- 
tion forces, from the corona region (the sphere with its 
radius r*) and the accretion disk, acting on the gas lo- 
cated at a location {p) in the field are assumed to be only 
in radial direction. The magnitude of the radiation force 
due to the corona is assumed to be a function of radius 
r, but that due to the accretion disk is assumed to be a 
function of r and the angle {9') between the disk angu- 
lar momentum Jd and the position vector r as shown 
in the figure. The point-source like approximation for 
the disk radiation pressure at P is valid when td 
In the following, we will describe our radiation hydrody- 
namics, our implementation of the continuum radiation 
sources (the corona and disk), the model parameters and 
assumptions. 

2.1. Hydrodynamics 

We employ 3-D hydro dynamical simulations of the out- 
flow from and accretion onto a central part of AGN, us- 
ing the ZEUS-MP code (c.f., Hayes et al. 2006) which 
is a massively parallel MPI-iniplenic!nted version of the 
ZEUS-3D code (c.f., Hardee & Clarke 1992; Clarke 1996). 
The ZEUS-MP is a Eulerian hydrodynamics code which 
uses the method of flnite differencing on a staggered 
mesh with a second-order-accurate, monotonic advection 
scheme (Hayes et al. 2006). To compute the structure 
and evolution of a flow irradiated by a strong contin- 
uum radiation of AGN, we solve the following set of HD 
equations: 



= - VP + pg + pg,^, 
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-PV -v + pC, 
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(3) 



where p, u, P and v are the mass density, energy 
density, pressure, and the velocity of gas respectively. 
Also, g is the gravitational force per unit mass. The 

Lagrangian/co- moving derivative is deflned as D/Dt = 
d/dt + v-'V. We have introduced two new components to 
the ZEUS-MP in order to treat the gas dynamics more 
appropriate for the gas flow in and around AGN. The 
first is the acceleration due to radiative force per unit 
mass (fifrad) ill equation (2), and the second is the the ef- 
fect of radiative cooling (and heating) simply as the net 
cooling rate (C) in erjuation (3). As this is our first 3-D 
simulations with this code, we consider a simplest case 
i.e., C = , but Qrad 7^ 0. We also use 7 = 1.01 in the 
equation of state P = (7 — 1) u where 7 is the adiabatic 
index. In the following, our implementation of g^^^^ will 
be described. 



2.1.1. Radiation Force 

To evaluate the radiative acceleration due to line ab- 
sorption/scattering, we follow the method in Proga et al. 
(2000) who applied the modified Castor, Abbott & Klein 
(CAK) approximation (Castor et al. 1975). Their model 
works under the assumption of the Sobolev approxi- 
mation (e.g., Sobolev 1957; Castor 1970; Lucy 1971); 
hence, the following conditions are assumed to be valid: 
(1) presence of large velocity gradient in the gas flow, 
and (2) the intrinsic line width is negligible compared 
to the Doppler broadening of a line. Following Proga 
et al. (2000), the radiative acceleration of a unit mass at 
a point r can be written as 



9rad= / + 

Jn 



ndD, 



(4) 



where I is the frequency-integrated continuum intensity 
in the direction n, and Q is the solid angle subtended 
by the source of continuum radiation. Also, (Te is the 
electron scattering cross section. The force multiplier Ai 
is a function of optical depth parameter r which is similar 
to the Sobolev optical depth (c.f. Rybicki & Hummer 
1978), and can be written as 



r = 



CTePVth 

\Q\ 



(5) 



where Q = dvi/ dl is the directional derivative of the ve- 
locity field in direction n, dl is the line element in the 
same direction, and Vth is the thermal velocity of the gas. 
Further equation (4) can be simplified greatly when the 
continuum radiation source is approximated as a point, 
i.e., when r ^ Vc where Vc is the radius of the radiation 
source. In our case, we consider the accretion disk which 
emits most of the radiation from the innermost part, be- 
tween r* and td in Figure 1; hence, the condition r ^ rc 
is satisfied. Using this approximation, the radiative ac- 
celeration gij.ad will be radial only, and be a function of 
radial position and polar angle (if the contribution from 
the disk luminosity is included), i.e. g^^d — 5rad {r,9)r. 
This simplification is very useful for our purposes as it 
reduces the computational time significantly hence it en- 
ables us to perform large-scale 3-D simulations. Unlike 
Proga et al. (2000) , we consider the case in which the ra- 
diative acceleration is dominated by the continuum pro- 
cess, i.e. = in equation (4) in this paper since we 
initially intend to investigate the basic characteristics of 
the impact of the disk precession that do not depend on 
the details of the radiation force model. The models with 
A4 7^ in equation (4) and C 7^ will be presented in a 
forthcoming paper. 

2.2. Continuum Radiation Source 

As mentioned earlier, we consider two different contin- 
uum radiation sources in our models: (1) the accretion 
disk, and (2) the central spherical corona. Since the ge- 
ometry of the central engine in AGN is not well under- 
stood, we assume that is consist of a spherically shaped 
corona with its radius and the innermost part of the 
accretion disk (c.f. Fig. 1). The disk is assumed to be 
flat, Keplerian, geometrically- thin and optically thick. 
The disk radiation is assumed to be dominated by the 
radiation from the disk radius between r* and td where 
r* = 3rs and r* < td <C n (c.f.. Fig. 1). Note that the 
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Fig. 1. — Basic model configuration. A super massive blackhole (BH) with its Schwarzschild radius rs is located at the center of the 
cartesian coordinate system (a;, y, z) where the y-axis is perpendicular to and into the page. The normal vector or the angular momentum 
{Jr>) of the accretion disk, spanning from its inner radius r* to its outer radius r^j, is misaligned with the z-axis by a tilt angle i.e., the 
accretion disk is tilted by from the x-y plane. The accretion disk hence its angular momentum Jd is assumed to precess around the 
z-axis with the precession period P. The 3-D hydrodynamic sinmlatioiis are performed in tlic splicrical coordinate system (r, 6, (f). The 
simulations are performed in the radial range between the inner boundary r\ and the outer boundary Tq. The radiation pressure from the 
central BH on a point p with its position vector r is in radial direction, and is a function of r. Whereas the radiation pressure from the 
accretion disk is assumed to be a function of r and 9' where the latter is the angle between Jjj and r (see Sees. 2.2 and 2.3 for details). 
The point-source like approximation for the disk radiation pressure at P is valid when <Sir\. Note that the figure is not to scale. 



exact size of ro docs not matter as long as it satisfies 
this condition in order for the point-source approxima- 
tion mentioned in §. 2.1 to be valid. 

In terms of the disk mass-accretion rate (Md), the 
mass of the BH (Mbh) and the Schwarzschild radius (rs), 
the total luminosity (L) of the system can be written as 



2r]GMBiiMr) 



(6) 
(7) 



where 77 is the rest mass conversion efficiency (e.g., 
Shakura & Sunyaev 1973). Following Proga (2007) and 
Proga et al. (2007), we simply assume the system es- 
sentially radiates only in the UV and the X-ray bands. 
The total luminosity of the system L is then the sum 
of the UV luminosity iuv and the X-ray luminosity Lx 
i.e., L = L\jY + Lx- Further, wc assume that the disk 
only radiates in the UV and the central corona in the 
X-ray. The ratio of the disk luminosity (^d) to the to- 
tal luminosity is parametrized as /d = Ly)/L, and that 
of the corona luminosity (L«) to the total luminosity as 
/* =L^IL. Consequently, /d + /* = 1- 



In the point-source approximation limit, the radiation 
flux from the central X-ray corona region can be written 
as 

^' = iS? 

where r is the radial distance from the center (by ne- 
glecting the source size). Here we neglect the geometri- 
cal obscuration of the corona emission by the accretion 
disk and vice versa. On the other hand, the disk radia- 
tion depends on the polar angle because of the source 
geometry. Again following Proga (2007) and Proga et al. 
(2007) (see also Proga ct al. 1998), the disk intensity /d 
is assumed to be radial and Id oc |cos0'|. This follows 
that the disk radiation flux at the distance r from the 
center can be written as 



= 2 |cos6i'| 



47rr2 



(9) 



where Q' is the angle between the disk normal and the 
position vector r (c.f.. Fig. 1). The leading term 2 in 
this expression comes from the normalization of the polar 
angle dependency. Finally by using eqs. (4), (8) and (9), 
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the radiative acceleration term in equation (3) can be 
written as 

9rad = ^{/*+2|cose'|/D}r. (10) 

2.3. Precessing Disk 

As we noted before, here we do not model the preces- 
sion of the accretion disk itself, but rather manually force 
the precession. We do not specify the cause of the preces- 
sion either. We simply assume that the disk precession 
exists, and investigate its consequence to the AGN en- 
vironment. The UV emitting portion of the disk spans 
from r* to ro (c.f.. Fig. 1). We assume that ru ^ n 
where the n is the inner radius of the computational do- 
main of the hydrodynamic simulations. This means that 
the disk itself is not in the computational domain. The 
effect of the processing disk is included as processing ra- 
diation field in the hydrodynamics of the gas (through 
eq. [2]). 

Wo assume that the disk is tilted from the x-y plane 
(in the cartesian coordinate system) by G as in Figure 1. 
Equivalently, the disk angular momentum Jd (assuming 
a flat uniform Keplerian disk) deviates from the z-axis by 
O. Further, we assume that Jd precesses around the z- 
axis with precession period P. With these assumptions, 
the components of the Jd in the cartesian coordinate 
system can be written as 

Jdx = Jd sin 6 cos (-^^ > (H) 

Jny = Jd sin 9 sin ^-^^ , (12) 

JD^ = Jr.cose (13) 

where t is the time measured from the beginning of hy- 
drodynamic simulations. Here we set Jd to be on the 
X — z plane (as shown in Fig. 1) at t = 0. By setting 
G = 0, the model reduced to an asymmetric case as in 
Proga (2007). To compute the radiative acceleration as 
expressed in equation (10), one requires the angle be- 
tween Jd and the position vector r at which the set of 
the HD equations (eqs. [1], [2] and [3]) are solve. This 
can be obtained simply by finding the inner product of 
Jd and r. 

2.4. Model Setup 

In all models presented here, the following ranges of 
the coordinates are adopted: ri < r < Tq, < < tt and 
<(}()< 27r where n = 500 r, and = 2.5 x 10^ r*. The 
radius of the central and spherical X-ray corona region r* 
coincides with the inner radius of the the accretion disk 
(Fig. 1). In our simulations, the polar and azimuthal 
angle ranges are divided into 128 and 64 zones, and are 
equally spaced. In the r direction, the gird is divided 
into 128 zones in which the zone size ratio is fixed at 
Arfe+i/Arfe = 1.04. 

For the initial conditions, the density and the tem- 
perature of gas are set uniformly i.e., p = po and 
T = To everywhere in the computational domain where 
Po = 1.0 X 10-21 gcm-3 and = 2 x lO"^ K through out 
this paper. The initial velocity of the gas is simply set 
to zero everywhere. 

At the inner and outer boundaries, we apply the out- 
flow (free-to-outflow) boundary conditions, in which the 



field values arc extrapolated beyond the boundaries us- 
ing the values of the ghost zones residing outside of nor- 
mal computational zones (see Stone & Norman 1992 for 
more details). At the outer boundary, all HD quantities 
(except the radial velocity) -atc fixed constant, to their 
initial values (e.g., T = To and p = po), during the the 
evolution of each model. The radial velocity components 
arc allowed to float. Proga (2007) applied these condi- 
tions to represent a steady flow condition at the outer 
boundary. They found that this technique leads to a 
solution that relaxes to a steady state in both spheri- 
cal and non-spherical accretion with an outflow (sec also 
Proga & Begelman 2003) . This imitates the condition in 
which a continuous supply of gas is available at the outer 
boundary. 

3. RESULTS 
3.1. Reference Values 

We consider four different cases which have different 
combinations of the disk tilt angle (G) and the disk pre- 
cession period (P), as summarized in Table 1. The fol- 
lowing parameters are common to all the models pre- 
sented here, and are exactly the same as in Proga (2007). 
We assume that the central BH is non-rotating and has 
mass A'/bh = 10*^ Mq. The size of the disk inner radius 
is assumed to be r* = 3rs = 8.8 x 10^^ cm (c.f. Sec. 2.4). 
The mass accretion rate (Ma) onto the central SMBH 
and the rest mass conversion efficiency (r/) arc assumed to 
be 1 X 10^^ g s^i and 0.0833, respectively. With these pa- 
rameters, the corresponding accretion luminosity of the 
system is L = 7.5 x lO'^^ergs-^ = 2 x lO^^L©. Equiv- 
alently, the system has the Eddington number F = 0.6 
where F = L/L^dd and Lsdd is the Eddington lumi- 
nosity of the Schwarzschild BH i.e., AwcGMBii/cre- The 
fractions of the luminosity in the UV (/uv) and that in 
the X-ray (/x) are fixed at 0.95 and 0.05 respectively, as 
in Proga (2007) (see their Run C). 

Important reference physical quantities relevant to our 
systems are as follows. The Compton radius, Rc = 
GA/bh/^ rrip/kTc, is 8 x 10"^* cm or equivalently 9 x 10'' 
where Tq, p and rup are the Compton temperature, the 
mean molecular weight of gas and the proton mass, re- 
spectively. Here we assume that the gas temperature at 
infinity is Too = Tc = 2 x 10^ K; hence, the correspond- 
ing speed of sound at infinity is Coo = {ikTc / prapY/"^ 
, x 10''cms-i where 7 is the adiabatic index. In this 
paper, 7 = 1.01 (almost isothermal) is adopted to im- 
itate a gas in Compton equilibrium with the radiation 
field. The corresponding Bondi radius (Bondi 1952) is 
Rb = GMbh/c^ = 4.8 X 10^* cm while its relation to the 
Compton radius is Rb = 'y~^Rc- The Bondi accretion 
rate (for the isothermal flow) is Mb = 3.3 x 10^^ gs-^ = 
O.52M0yr^i. The corresponding free-fall time (iff) of 
gas from the Bondi radius to the inner boundary is 
2.1 X 10" sec = 7.0 x lO^yr which is about 2.3 times 
smaller than the precession period used for Models II and 
HI, and about 23 times smaller than that of Model IV 
(c.f. Tab. 1). 

3.2. Comparison of axisymmetric models in 2-D and 

3-D 

Before we proceed to the main precession disk models, 
we briefly compare our axisymmetric model (Model I) 
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TABLE 1 
Model Summary 



Model (°) (yr) (lO^^gs-^) (lO^^gs-^) (lO^^gs-^) (kms'^) (lO^"ergs-i) (lO^°ergs-^) (jo) 

I oo -2.2 -0.6 1.6 1500 3.0 260 

II 5 1.6 X 10* -2.3 -0.6 1.7 640 1.4 290 0.2 

III 15 1.6 X 10"* -2.2 -0.4 1.8 620 1.1 310 0.2 

IV 5 1.6 X 10^ -2.2 -0.6 1.6 1900 3.0 260 0.05 



The model output values are averaged over the last 2 x 10^^ s of the hydrodynamic simulations. 



Note. — 

with the axisymmetric models presented earlier by Proga 
(2007) who used very similar model parameters as in our 
Model I. The main differences here are in the treatment 
of the radiation force and that in the radiative heat- 
ing/cooling. As mentioned earlier, we set the force mul- 
tiplier M. = Q {in eq. [4]) and the net cooling rate C = 
(in eq. [3]) while Proga (2007) used non-zero values of 
those two terms. In our Model I, the adiabatic index is 
set to 7 = 1.01 (essentially isothermal), but their models 
use 7 = 5/3. However, Proga (2007) found that their 
Run A is nearly isothermal despite 7 = 5/3 was used 
(see their Fig. 1). Another important difference is the 
numerical codes used. Proga (2007) used the ZEUS-2D 
code (Stone & Norman 1992). 

Overall geometry of the flow in Model I (Figs. 2 and 
3) is similar to those in Proga (2007). The matter ac- 
cretes onto the central BH near the equatorial plane, and 
strong outflows occur in polar direction. The collimation 
of our model is relatively weak compared to their Run C 
which uses exactly the same disk and corona luminosities 
as in our Model I. The wider bipolar outflow pattern seen 
here resembles that of their Run A which has the high- 
est X-ray heating. The difference and the resemblance 
seen here arc caused by the following two key factors: 
(1) nearly isothermal equation of state and (2) no radia- 
tive cooling (C = 0) in our model . These condition keep 
the temperature of gas warm everywhere in the compu- 
tational domain, and the temperature is essentially that 
set at the outer boundary (Too = Tb = 2 x lO'^ K). This 
will result in a very similar situation as in Run A of 
Proga (2007) in which the gas temperature is also rela- 
tively high because of the high X-ray heating and cooling. 
The high temperature hence the ionization state of the 
gas makes the line force in their model very inefficient, 
resulting in the situation in which the gas is almost en- 
tirely driven by the continuum process (electron scatter- 
ing) and thermal effects just as in our Model I. 

Although not shown here, we have also checked the 
internal consistency of the ZEUS-MP (3-D) code by run- 
ning the axisymmetric models (Model I) in both 2-D and 
3-D modes. We find that the results from the both runs 
agree with each other in all aspects e.g., inflow and out- 
flow geometry, density distribution, velocity, mass accre- 
tion and outflow rates. 

3.3. Dependency on the disk tilt angle © 

We now examine the model dependency on disk tilt 
angle (8) while keeping all other parameters flxed. The 
results from Models I, II and III (c.f.. Tab 1), which 
use 6 = 0°, 5° and 15° respectively, are compared for 
this purpose. Note that the observations suggest that a 
large fraction of AGN have rather small i.e., <'~ 15° 



(e.g., Lu & Zhou 2005). 

Figure 2 shows the slices of the density and velocity 
flelds (on the y = plane) from snapshots of our four 
simulations. The snapshots are chosen at the time when 
the models reached a (semi-)steady state for Models I 
and II. As we will sec later, the flow never reaches steady 
state in Model III; therefore, we chose the snapshot of 
the model at the time when the flow pattern is a typ- 
ical of a whole simulation time sequence. While accre- 
tion occurs mainly on the equatorial plane [z — 0) for 
Model I, it occurs in a inclined plane with a pitch angle 
(a angle between the equatorial plane and the accretion 
plane) similar to the disk inclination angle, for Models II 
and III. In the precessing disk models (II and III), the 
deviation from the axisymmetric is clearly seen in both 
density distribution of gas and the shapes of the Mach 
number contours. Corresponding 3-D density and Mach 
number contour surfaces of these models are also shown 
in Figs. 3 and 4. The morphology of the density distri- 
bution seen in Figure 2 resemble that of the Z-shaped 
(for Model II) and the S-shaped (for Model III) radio 
jets (e.g. Condon & MitcheU 1984; Hunstead et al. 1984; 
Tremblay et al. 2006) although in different scales. Ob- 
viously, the difference between the Z- and S- shapes are 
simply due to the difference in the viewing angles. Unlike 
the MHD precessing jet models, the bending structures 
of the density distributions seen here are shaped by the 
geometry of the sonic surfaces. When accreting mate- 
rial from the outer boundary encounters the relatively 
low density but high-speed outflowing gas launched by 
the radiation force from the inner part, the gas becomes 
compressed, and forms higher density regions. The flows 
in the bending density structure itself are rather com- 
plex (especially in Model III), but the direction of the 
flow becomes outward (in radial direction) as they ap- 
proach the sonic surface (excluding the one shaped like 
a disk formed by the accreting gas in the inner region). 
Relatively large curvatures of the flows seen in both den- 
sity and the Mach number contours of Models II and III 
can be also understood from the fact the precession pe- 
riod used in these models [P = 16000 yr) is comparable 
to the gas free fall time (tg = 7000 yr, c.f., § 3.1). The 
curvatures or the "twists" of the weakly coUimated bipo- 
lar flows can be clearly seen in the 3-D representation of 
these models in Figs. 3 and 4. 

We compute the mass fluxes as a function of radius 
for a quantitative comparison of the characteristics of 
the flows in the models. Following Proga (2007), the 
net mass flux (Mnet), the inflow mass flux (Mi„) and the 
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Fig. 2. — Comparison of the density and velocity fields on the x-z plane from Models I (upper-left) II (upper-right), III (lower-left) and 
IV (lower-right). The density maps shown in the background are given in logarithmic scale (base 10) and in cgs units. The contours of the 
Mach number are overlaid along with the arrows which indicate the directions of the velocity on the x-z plane. The units of both x and z 
axes are in pc. The time slices of each models are chosen such that the density and velocity fields are representative of each model. 



outflow mass flux (-Mout) can be computed from 



M{r) 



V ■ da 



(14) 
(15) 



where Vr is the radial component of velocity v. In the 
equation above, M = -Mnet if all Vr are included. Sim- 
ilarly, M = Min for Vr < and M = Mout for Vr > 0. 
Also, da — f sin 9 dd dip and dil = sin 6 d9 d(j). Simi- 
larly we further define the outflow power in the form of 
kinetic energy (Pk) and that in the thermal energy (Pth) 



as functions of radius i.e., 



and 



Pk{r) 



^th (r) 



pv; 



dn 



.dn. 



(16) 



(17) 



where u,. > 0. 

The resulting mass fluxes and the outflow powers of 
the models are summarized in Figure 5. In all cases, 
the mass inflow flux exceeds the mass outflow rate at all 




Fig. 3. — Two-level iso-density (left panels) and the corresponding sonic (right panels) surfaces in 3-D for Models I (upper panels) and II 
(lower panels). The density levels used here are log p = —20.5 (blue) and —21 (red) where p is in gcm~'^. The time slices of the simulation 
data used here are as in Fig. 2. The sizes of the plotting boxes are 14.2 pc in all directions {x, y and z). 



radii. For Models I and II, the net mass fluxes (Mnot) 
are almost constant at all radii, indicating that the flows 
in these models are steady. Despite the presence of the 
disk precession in Model II, the flow becomes steady. 
The density distribution and the velocity field become al- 
most constant in the coordinate system co-rotating with 
the disk precession period. On the other hand, Afnct for 
Model III does not remain constant as r becomes larger 
(r > 10^* cm) because of the unsteady nature of the flow 
(c.f.. Figs. 2 and 4). As the disk tilt angle Q increases, 
the direction of the outflows, which are normally in po- 
lar directions (±z directions) with an absence of the disk 
tilt, moves toward the equatorial plane (the x-z plane) 



where the flow is predominantly inward. This opposite 
flows makes it harder for the outflowing gas to reach the 
outer boundary. Further, since the disk is processing, the 
direction of the outflow is constantly changing. This re- 
sults in continuous collisions between the inflowing and 
outflowing gas especially for a larger Q model. The net 
mass fluxes at the inner boundary M^ct (n) are —0.6, 
-0.6 and -0.4 x lO^^ gs'^ (or equivalently -0.10, -0.10 
and —0.06 M© yr~^) for Models I, II and III respectively 
(Tab. 1), indicating the net mass flux inward (negative 
signs indicate inflow) decreases slightly, but not signifi- 
cantly as the disk tilt angle 8 increases. 
The ratios of the total mass outflow flux to the total 






Fig. 4. — As in Fig. 3, but for Models III and IV. 



mass inflow at the outer boundary (/i — Mout/Mj 

are 0.73, 0.74, 0.82 for Models I, II, and III (see also 
Tab. 1). These values indicates that the high efficiency 
of the outflow by the radiation pressure even for a modest 
Eddington number used here i.e., T = 0.6. This conver- 
sion efficiency /i (from the outflow to inflow) is about 
the same for Models I and II, but it slightly (~ 12%) 
increases for Model III which has the highest disk tilt 
angle. Overall characteristics of the mass-flux curves as 
a function radius for Models I and II are also very sim- 
ilar to each other. The curves for Model III are also 
similar to those of Models I and II; however, they differ 
in the outer radii (r <~ lO^^cm), mainly because of the 



unsteady nature of the flow in this model. 

The maximum speed of the outflow in the radial 
direction w™'^'^ (ro) decreases as 9 increases (Tab. 1). 
The reduction in the speed is very signiflcant {Avr = 
— 860kms~^) as 9 increase from 0° to 5° while the 
change is relatively small {Avr = — 20kms^^) as 9 
changes from 5° to 15°. 

Figure 5 also shows the outflow powers {Pk and Pth) of 
the models as a function of radius, as defined in eqs. (16) 
and (17). As for the mass flux curves in the same fig- 
ure, the dependency of the energy flux curves on radius 
for Models I, II and III are very similar to each oth- 
ers. A small but noticeable deviations of the curves for 
Model III from those for Models I and II are seen at 




Fig. 5. — Comparison of mass flux and energy flux as a function of radius for Models I (upper- left) , II (upper-right), III (lower-left) 
and IV (lower-right). The plot for each model is subdivided into two panels: top (mass flux) and bottom (energy flux). In the mass flux 

plots, the inflow (dash: Min), outflow (solid: Mo) and net (dot: A^not) mass fluxes, as defined in equation (15), are separately plotted as a 
function of radial distance from the center. The absolute values of Min and Mnet are plotted here since they are negative at all radii. In 
the energy flux plots, the kinetic energy (solid) and the thermal energy (dot), deflned in eqs. (16) and (17), are shown. Note that the times 
slices of the model simulations used here to computed the fluxes are same as those in Figure 2. 



the large radii (r > 3 x 10^^ cm). The figure shows that 
in all three models, the outflow power is dominated by 
thermal process (Pth « 10 — 100 at all radii). This 
can be explained by the high temperatures of the gas 
[T Ki Tc = 2.0 X 10^ K ) in the computational domains 
caused by the (almost) isothermal equation of state and 
the temperature (2.0 x 10^ K) fixed at the outer bound- 
ary. The kinetic powers or the radiation forces are not as 
significant as the pressure gradient force in these models; 
however, their importance cannot be ignored since they 
"shape" the geometry of the outflow as they strongly de- 
pend on the polar angle position of a point in the com- 



putational field. We also note that as © increases, the 

kinetic power at the outer boundary Pk {to) decreases sig- 
nificantly e.g. Pk {tq) of Model III is three times smaller 
than that of the axisymmetric model, Model I (see Ta- 
ble 1). 

Next we examine the degree of non-axisymmetry in 

Model I, II and III, and seek for any obvious dependency 
on G. For this purpose, we compute the center of mass 
(CM) of the gas on the planes perpendicular to the z- 
axis (as a function of z) i.e., Xc {z) and j/c {z) which are 
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defined as 

and 

Vc {z) 

where 
m (z) 



f-r 11a I^^oo x5{z- z')p (x, y, z') dz' dx dy 



m (z) 



(18) 



I-r I-a I^oo VH^- z')p{x, y, z') dz' dx dy 



m (z) 



(19) 



(5 (z — z') p (a;, y, z') dz' dx dy 



(20) 

1/2 

and a {y, z) — (r^ — — y^) . The results are shown 
in Figure 6. As expected, the CM position remains con- 
stant and on the z-axis {xc (z) = and j/c {z) = 0) for 
Model I, as this is an axisymmetric model. For both 
Models II and III, the maximum amount of deviations 
for each component of the CM (|a;c|and \yc\) is about 
0.3 pc which is relatively small compared to the outer 
boundary radius (ro = 7.1 pc). The Xc and yc curves 
are anti-symmetric about the z = position since our 
model accretion disk hence the radiation force is sym- 
metric about the origin of the coordinate system. The 
plot also shows that the positions of the maxima and 
minima in the Xc and yc curves do not coincide, but they 
are rather shifted in both -\-z and — z directions. This 
clearly demonstrates a helical or twisting nature of the 
flows, as one can also simply see it in the 3-D density 
and Mach number contour plots in Figs. 3 and 4. 

To summarize, as the tilt angle of the disk precession 
Q increases, reductions of the maximum outflow veloc- 
ity (Vr) and the kinetic outflow power (P^) at the outer 
boundary Tq occur, as a consequence of the stronger in- 
teractions between the outflowing and inflowing gas of as 
increases. The net mass inflow flux (Mnot) at the in- 
ner boundary does not strongly depend on Q. The ther- 
mal outflow energy power dominates the kinetic outflow 
power in our models here because of the high tempera- 
ture of set at the outer boundary and because the gas 
is (almost) isothermal. The flows of Models II and III 
show helical structures; however, the radius of the he- 
lices (base on the CM positions along the z-axis) does 
not change greatly as Q increases from 5° to 15°. 

3.4. Dependency on the disk precession period P 

We now examine the dependency of the model on the 
disk precession period (P). We vary the value of P while 
flxing the disk tilt angle to = 5°. For this purposed, we 
compare Models I, II and IV as summarized in Table. 1. 
The precession periods P are oo, 1.6 x 10** and 1.6 x 10^ yr 
respectively for Models I, II and IV. In the units of the 
free-fall time (ig — 7.0 x 10'^ yr ) from the Bondi radius 
(§ 3.1), they are oo, 2.3 and 23 respectively. Note that the 
observations suggest that typical values of jet precession 
period are P = lO-^-lO^ yr (c.f., Tab.l in Lu & Zhou 
2005). 

Figure 2 shows that similarities between Model IV and 
Model I in their morphology of the density distribution 
and Mach number surfaces. At a given time, the flow in 
Model IV is almost axisymmetric, and the symmetry axis 
is tilted also by = 5° from the z-axis. This is caused by 
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Fig. 6. — Positions of the center of mass of the gas on the planes 
perpendicular to the 2-axis for Models I, II, III and IV from the 
top bottom. The locations of the center of mass, Xc (sold) and j/c 
(dash) as defined in eqs. (18) and (19), are computed as a function 
of z. All models except Model I show a clear sign of deviations 
from axisymmetry; however the displacement of the center of mass 
remains fairly small (|xc| <~ 0.3pc and Ij/d <~ 0.3pc), at all 
z locations, compared to the size of outermost radius To (7.1 pc) 
of the computational domain. The patterns in Xc and yc curves 
for Models 11 and III indicate that the flow density structures are 
helical. 



the relatively long precession period for Model IV com- 
pared to the dynamical time scale or the gas free-fall 
time scale iff. The curvature or helical motion of the 
gas is not significant, and it does not greatly affect the 
overall morphology of the flow, except for the outermost 
part of the flow where the flow is slightly turbulent due 
to the shear of the slowly precessing flow and the outer 
boundary. This can be clearly seen in the 3-D plots Fig- 
ure 4. As the precession period P becomes shorter and 
comparable to iff (as in Model II), the flow shows more 
curvature and the helical structures. 

The mass flux curves (c.f., eq. [15]) for Models IV in 
Figure 5 also show that nature of the flows between Mod- 
els I and IV are very similar to each other. Overall 
characteristics of the curves are also similar to that of 
Model II. In fact, the net mass flux Mnet, the inflow mass 
flux Mill and the outflow mass flux Mout at the outer 
boundary of Model IV are identical to those of Model I 
(see Table 1). Also note that Models I, II and IV all have 
same Min value at the inner boundary, i.e., the mass in- 
flow rate across the inner boundary in insensitive to the 
change in the precession period for = 5°. 

The similarity between and Models I and IV can be also 
seen in the outflow powers, and Pth. Figure 5 shows 
Pfc and Pth as a function of radius for Model IV are al- 
most identical to those of Model I. The Pfe and Pth values 
at outer boundary are indeed identical (Table 1). A slight 
increase in the maximum outflow velocity at the outer 
boundary v^^^ija) is observed for Model IV, compared 
to Model I. The kinetic outflow power Pfc and the max- 
imum outflow velocity at the outer boundary v^'^^ira) 
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decreases as the precession period become comparable 
to iff. 

The CM positions Xc and Uc as a function of z (see 
§3.3) for Model IV is shown in Figure 6. Compared to 
Model II, the maximum displacement of the CM is about 
40 times smaller in Model IV i.e. |a;c| <^ 0.075 pc and 
\yc\ <^ 0.075 pc. The Xc curve for Model IV shows a 
rather complex pattern compared to that in Model II. 
This and the visual inspection of the density and the 
Mach number contour surfaces in Figure 4, indicates that 
the bipolar outflow flows are slightly twisted, but does 
not have clear helical structure. 

3.5. Time Evolution of Mass Accretion/ Outflow Rates 

and Angular Mom,entum 

Next, wc examine the variability or steadiness of the 
flows in each model by monitoring the mass fluxes at 
the outer boundary as in equation (15) and the angular 
momentum of the system as a function of time. For the 
latter, we compute the density-weighted mean speciflc 
angular momentum of the systems deflned as: 

_ jyp{rxv) dV 

'o- UpdV ^''^ 

where the denominator is simply the total mass of the 
gas in the computational domain. Note that the radi- 
ation force (eq. [10]) and the gravitational force are in 
radial direction only. Consequently, they do not exert 
torque onto the system; hence, they do not contribute 
to the change in the angular momentum of the system 
directly. The system can gain the angular momentum 
in the following way. In our models, the strength of the 
disk radiation field depends on the angle measured from 
the disk normal (c.f., eq. [9]). This causes gas pressure 
gradients in azimuthal direction, and contributes to the 
angular momentum of gas locally, forming vorticity. The 
precession of radiation field hence the precessing outfiow 
will cause the gas with preferred sign of vorticity to es- 
cape from the outer boundary, resulting a change in the 
net angular momentum of the gas in the computational 
domain. 

Figure 7 shows Mjn (r,,), Mout (to) and jp for Mod- 
els I, II, III and IV as a function of time. For Mod- 
els I, II and IV, both the mass fluxes and jp reach to 
asymptotic values by t « 7 x 10^^ s. Small oscillations 
of jp around the asymptotic values are seen for Mod- 
els II and IV. On the other hand, the mass fluxes of 
Model III has much larger amplitudes of the oscillations 
around an asymptotic value. By visual inspections, their 
oscillations do not seem have a clear periodicity associ- 
ated with them. We performed the Lomb-Scargle pe- 
riodgram analysis (e.g. Home & Baliunas 1986; Press 
et al. 1992) on the Min {ro), Mout {fo) and jp curves for 
Model III. Only Mout (?'o) shows a relatively strong signal 
at Pls = 1-36 X 10"^^ s which is about 2.7 times longer 
than the precession period of Model III. On the other 
hand A'/jn (ro) and jp curves do not have any obvious pe- 
riod associated with them, but they arc rather stochastic. 

As mentioned in § 3.3, as the disk tilt angle 6 increases 
the direction of the outflows, which normally exists in 
polar directions with an absence of the disk tilt, moves 
toward the equatorial plane (the x-z plane) where the 
flow is predominantly inward. In addition, the precession 



of the disk causes the direction of the outflow to change 
constantly; hence, causing constant creation of the shock 
between the inflowing and the outflowing gas. This leads 
into a very unstable flow of the gas at all time for a model 
with a larger 9 e.g.. Model III with 9 = 15°. The flow, 
of course, can be stabilized if the precession period is 
increased to a value much larger than the free-fall time 
iff. 

The amount of the (density-weighted) mean speciflc 
angular momentum deposited to the gas by the precess- 
ing disk (measured by jp) is largest in Models II and III 
(see Table 1), and that in Model IV is about 4 times 
less than those of Models II and III. For all models, a 
time-averaged value (by using the last 2 x 10"'^^ s of the 
simulation) oi jp is used in Table 1. It seems that the 
faster the disk precesses, the lager the amount of angular 
momentum transferred to the environment; however, this 
trend does not continues, as we increase the disk preces- 
sion speed even faster. Although not shown here, a model 
with exactly the same set of parameters as in Model II, 
but with P = 1600 yr (10 times faster rotation), showed 
that the the value of jp decrease to ^ 0.01, which is even 
smaller than Model IV (with P = 160000 yr). This indi- 
cates that the amount of angular momentum deposited 
to the gas depends on how close the precession period to 
the dynamical time scale of the flow. 

In principle, it is possible to model the change in the 
angular momentum of the accretion disk itself through 
the transfer of angular momentum from the environment, 
we ignored this effect for simplicity (and this is also our 
model limitation). To model the interaction of the disk 
angular momentum and the angular momentum of the 
surrounding gas properly, we need to model the dynam- 
ics of the gas in the accretion disk itself as well as the 
dynamic of the gas which is much larger scale as in our 
models here. This is computationally challenging with 
our current code since we have to resolve the length scale 
of the innermost part of the accretion disk (~ 10^^ pc) 
to the large scale outflow/inflow gas (~ 10 pc). 

4. CONCLUSIONS 

We have studied the dynamics of the gas under the in- 
fluences of the gravity of SMBH and the radiation force 
from the luminous accretion disk around the SMBH. The 
rotational axis of the disk was assumed to be tilted with 
respect to the symmetry axis with a given angle 9 and 
a precession period P (c.f., Figure 1). We have investi- 
gated the dependency of the flow morphology, mass ac- 
cretion/outflow rates, angular momentum of the flows 
for different combinations of 9 and P. This is a natural 
extension of similar but more comprehensive 2-D radia- 
tion hydrodynamics models of AGN outflow models by 
Proga et al. (2000), Proga (2007). As this is our first 
attempt for modeling such gas dynamics in full 3-D, we 
have used a reduced set of physical models described in 
Proga (2007) i.e., the radiation force due to line and dust 
scattering/absorption, and the radiative cooling/heating 
are omitted. In the following, we summarize our main 
findings through this investigation. 

(1) Our assumption of the adiabatic index (7 = 1.01) 
keeps the mean temperature of the gas in the compu- 
tational domain relatively high (~ 2 x lO^'K) which is 
essentially determined by the outer boundary condition. 
For our axisymmetric model (Model I: Figs. 2 and 3), 
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Fig. 7. — The mass flow rates across the outer boundary and the density-weighted mean specific angular momentum jp (c.f. eqs. [15] 
and [21]) plotted as a function of time for Models I (upper-left), II (upper-right), III (lower- left) , and IV (lower-right). For each model, the 
plot is subdivided into two panels: mass-inflow/outflow rates (top panels) and jp (bottom panels). In the top panels, the mass-inflow rate 
at the outer boundary (solid), and the mass-outflow rate at the outer boundary (dash) are shown separately. The values of jp are in the 
units of JO which is deflned as the speciflc angular momentum of the gas in Kcplerian orbit at the inner boundary {r = rin). Note that the 
precession period used here are oo, 5.0 x 10" s, 5.0 x 10" s and 5.0 x lO^^ s for Models I, II, III and IV respectively 



this results in the flow morphology very similar to the 
model with a relatively high X-ray heating (see Run A 
in Proga 2007) in which the line force is inefficient be- 
cause of the high gas temperature and hence the high 
ionization state of the gas. 

(2) Although in different scales, we were able to re- 
produced the Z- or S- shaped density morphology of the 
gas outflows (Fig. 2) which arc often seen in the radio 
observations of AGN (e.g. Florido ct al. 1990; Hutchings 
ct al. 1988). The bending structure seen here arc shaped 
by the shape of the sonic surfaces. When accreting ma- 
terial from the outer boundary encounters the relatively 
low density but high speed outflowing gas launched by 



the radiation force from the inner part, the gas becomes 
compressed, and forms higher density regions. 

(3) As the tilt angle of the disk precession increases, 
the reduction of the maximum outflow velocity [vj.) and 
the kinetic outflow power Pk at the outer boundary 
decrease as a consequence of the stronger interactions 
between the outflowing and inflowing gas (Tab. 1). The 
net mass inflow rate (Afnct) at the inner boundary does 
not change signiflcantly with increasing Q. 

(4) A relatively high efliciency of the outflow {^, = 
Mont/ Mill) by the radiation pressure were observed in 
our models (70-80 %; see also Tab. 1) for a Eddington 
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number (F) of 0.6 here. The conversion efficiency ji (from 
the outflow to inflow) is about the same for Models I and 
II, but it slightly (~ 12%) increased for Model III which 
has the highest disk tilt angle. 

(5) The thermal outflow energy power dominates the 
kinetic outflow power (Fig. 5) in the models presented 
here because of the high temperature of the flow (as men- 
tioned above). 

(6) The flows of Models II and III show helical struc- 
tures (c.f., Figs. 3 and 4); however, the radius of the 
helices does not change as 9 increases from 5° to 15°, 
based on the locations of the center of mass (Figure 6) 
of the planes perpendicular to the symmetry axis (the 
z-axis in Fig. 1). We leave for a future investigation to 
test whether these trends continue as becomes larger 
than 15°. 

(7) The characteristics of the flows are closely rclatcid 
to a combination of P and G, but not to P and O indi- 
vidually. Even with a relatively large disk tilt angle ©, 
if the precession period is much larger than the dynam- 
ical time scale of a system, the flow geometry obviously 
becomes almost axisymmetric (c.f. Model IV in Figs. 2 
and 4). 

(8) The gas dynamics of a model with a relative large 
disk tilt angle (6 = 15°) with a precession period 
comparable to the gas free-fall time (iff) of the system 
(e.g., Model III) does not reach a steady state because 
the outflows driven by the luminous accretion disk con- 
stantly collides with the inflowing/accreting gas as the 



disk processes hence as the outflow direction changes. 

(9) The amount of the density-weight mean specific 
angular momentum {jp) deposited by the precessing disk 
is largest for Models II and III (Tab. 1 and Fig. 7) which 
have the precession period comparable to tg. 

The models represented here are mainly for ex- 
ploratory purpose - to examine the basic model de- 
pendencies on Q and P - with a relatively simple set 
of physics but in full 3-D. In the follow-up paper, we 
will improve our model by including the physics omit- 
ted here (the line scattering/absorption, dust scatter- 
ing/absorption, and the radiative cooling/heating) as in 
the 2-D models of e.g., Proga (2007). 
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